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We present a novel approach to calculate the effective exchange interaction 
parameters based on the realistic electronic structure of correlated magnetic 
crystals in local approach with the frequency dependent self energy. The 
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analog of "local force theorem" in the density functional theory is proven for 
highly correlated systems. The expressions for effective exchange parame- 
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Q\ ' netic iron, with the local correlation effects from the numerically exact QMC- 
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The first-principle calculations of magnetic excitation spectrum for ferromag- 



scheme is presented. 
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I. INTRODUCTION 



The calculation of thermodynamic properties and excitation spectra of different magnetic 
materials is one of the most important motivations of the microscopic theory of magnetism. 
The main approach for such type of investigations is the local spin density functional (LSDF) 
scheme M. However, this method has some serious shortcomings when applying to transition 
metal and rare-earth magnetic materials. The main defect is the absence of the "Hubbard" 
type correlations which are most important for real magnets (see, e.g., recent reviews Pf-0]]). 
This leads to, generally speaking, incorrect description of electronic structure for such im- 
portant groups of magnetic materials as rare-earth metals and their compounds, metal-oxide 
compounds (including "classical" Mott insulators such as NiO and MnO as well as high-T c 
cuprates) and even for the iron-group metals ||,|4]|| . At the same time, the experience with 
the Hubbard model shows that the description of electronic structure and magnetic prop- 
erties of highly correlated materials are closely connected. Recently we propose a rather 
general scheme (so-called "LDA++ approach") for first-principle calculations of the elec- 
tronic structure with the local correlation effects being included 0. In this technique the 
full matrix of on-site Coulomb repulsion for the correlated states is taken into account in 
the local approximation for the electron self- energy. In such a way, we could provide a 
rather reasonable description of the electronic structure for different correlated systems such 
as Fe, NiO, and TmSe. It will be very useful to develop this approach for the description of 
different magnetic characteristics. 

The most rigorous way to consider properties of magnetic excitations is the calculation 
of frequency- dependent magnetic susceptibilities [|J . However, for many important cases we 
can restrict ourselves to more simple problem of the calculation of static response functions. 
More explicitly, one can consider the variations of total energy E to t (or thermodynamic 
potential Q) with respect to magnetic moments rotations. This approach results in the 
magnetic interactions of different types: the variation of total energy of a ferromagnet over 
the rotation of all spins at the same angle determines the magnetic anisotropy energy, while 
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the variation of E to t over the relative rotations of spins on two sites gives the parameters of 
pairwise exchange interactions, etc. This approach was proposed earlier in the the framework 
of the LSDF scheme [J]. It is sufficient for the calculation of "phenomenological" exchange 
parameters which are important for the consideration of domain wall widths and other 
"micromagnetic" properties. In the adiabatic approximation when the energies of magnetic 
excitations are small in comparison with typical electronic energies this is also sufficient for 
the calculation of the spin-wave spectrum. In the mean-field approximation these "exchange 
parameters" can be used for the estimation of Curie or Neel temperature. 

In this work we derive general expressions for the parameters of magnetic interactions 
in LDA++ approach and calculate the exchange parameters for ferromagnetic iron. It is a 
first attempt to investigate magnetic interactions, taking into account correlation effects in 
the electronic structure for real materials. 



II. GENERAL FORMALISM 

A. Local force theorem in LDA-| — |- approach 

An important trick for the definition of exchange interactions in the LSDF approach is 
the use of so called "local force theorem". This reduces the calculation of the total energy 
change to the variations of one-particle density of states |||7J . First of all, let us prove 
the analog of local force theorem in the LDA++ approach. In contrast with the standard 
density functional theory, it deals with the real dynamical quasiparticles defined via Green 
functions for the correlated electrons rather than with Kohn-Sham "quasiparticles" which 
are, strictly speaking, only auxiliary states to calculate the total energy. Therefore, instead 
of the working with the thermodynamic potential Q as a density functional we have to start 
from its general expression in terms of an exact Green function || 

n = Q sp — Vide (1) 
n sp = -Tr {in [S - G l ] } 



Q dc = TrEG - $ 

where G, Go and E are an exact Green function, its bare value and self-energy, correspond- 
ingly; $ is the Luttinger generating functional (sum of the all connected skeleton diagrams 
without free legs), Tr = Tr uiLa is the sum over Matsubara frequencies Tr u ... = TJ2---, 
uj = 7rT (2n + 1) , n = 0, ±1, T is the temperature, and iLa are site numbers (i), orbital 
quantum numbers (L = l,m) and spin projections a , correspondingly. We have to keep in 
mind also Dyson equation 

G- 1 = Gq 1 - E (2) 

and the variational identity 

5$ = TrTjdG. (3) 

We represent the expression (|I|) as a difference of "single particle" (sp) and "double counted" 
(dc) terms as it is usual in the density functional theory. When neglecting the quasiparticle 
damping, Q sp will be nothing but the thermodynamic potential of "free" fermions but with 
exact quasiparticle energies. Suppose we change the external potential, for example, by 
small spin rotations. Then the variation of the thermodynamic potential can be written as 

5Q = S*il sp + 5iil sp — 5ildc (4) 

where 5* is the variation without taking into account the change of the "self-consistent 
potential" (i.e. self energy) and 8% is the variation due to this change of E. Taking into 
account Eq. (|3|) it can be easily shown (cf. ||) that 

= 5n dc = TrG5Tj (5) 

and hence 



= 6*tt sp = -S*Tr In 



E - G, 1 } (6) 



which is an analog of the "local force theorem" in the density functional theory J7J. In 
the LSDF scheme all the computational results expressed in terms of the retarded Green 



function and not in the Matsubara one. The relations of "real" and "complex" Green- 
function formulae are based on the identity 

1 °° 

T ruJ F(iuj) = / dzf(z)ImF(z + iO) (7) 

7T J 



-00 



where f(z) = [exp z/T + 1] 1 is the Fermi function, F(z) is a function regular in all the 

complex plane except real axis. Therefore Eq.(^) takes the following form 

00 

5Q=- f dzf{z)ImTr iLa \nG- 1 {z + iO) (8) 

71 J 

—00 

which is the starting point for the calculations of magnetic interactions in LSDF approach 
0. However note, that in the case of frequency-dependent self-energy (LDA++ approach) 
it is more convenient to work with the Matsubara Green functions [[J. 

B. Effective spin Hamiltonian 

Further considerations are similar to the corresponding ones in LSDF approach. The 
most suitable way based on the sum rule is proposed in ||. 

In the LDA++ scheme, the self energy is local, i.e. is diagonal in site indices. Let us 
write the spin-matrix structure of the self energy and Green function in the following form 

= + EJo- (9) 
where E^'^ = | Csj ± E|V S| = Efej, with ej being the unit vector in the direction 



of effective spin-dependent potential on site i, a — {a X) a y) a z ) are Pauli matrices, G?- = 
^Tr a (Gij) and G*- = \Tr a {Gija). We assume that the bare Green function G° does not 
depend on spin directions and all the spin-dependent terms including the Hartree-Fock terms 
are incorporated in the self energy. Spin excitations with low energies are connected with 
the rotations of vectors 

5et = Sifi x ej (10) 



According to the "local force theorem" @ the corresponding variation of the thermodynamic 
potential can be written as 

5tt = 5*tt sp = ViSifi (11) 

where the torque is equal to 

V, = 2Tr UJ L [E{ x GJJ (12) 

Further we have to use an important sum rule for the Green function which is the consequence 
of Dyson equation. Using Eq. (0) one has 

G = G C -T+G S ■a = GG- 1 G= (13) 
G ((Gq 1 - S c ) • I - E s • a) (G c ■ I + G s ■ a) . 



Separating the spin-dependent and the spin-independent parts in Eq. ([X3|) we have the 
following sum rules for G c 



G c = G c (Go 1 - S C J G c - G C E S G S = G c (G 1 - £ c ) G c - G S E S G C (14) 

and similarly for G s 

G s = - (G S E S ) G s + G C E S G C + iG c (E s x G s ) = (15) 
= -G s (E S G S ) + G C E S G C + i (G s x E s ) G c 

Then for diagonal elements of the Green function one obtains 

Gu = - E [(g^) g;, - g;^g), - iG% (e- x gj<)] (ie) 
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Substituting ( |T5| ) into (|T2|) we have the following expression for the torque 



V, = 2Tr wi [SJ x G«] (17) 

= -2 £ •/ /•,, {e; x (g;e;) g% -e;x c%e;g& - ,e-; x G' ; (2:5 x G%) } 

If we represent the total thermodynamic potential of spin rotations or the effective Hamil- 
tonian in the form 
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n spin = -E {(g&ej) (g; 4 sj) - s|g°.e^ - ; (s| x g^e*) g;} (is) 



one can show by direct calculations that 

~5Q 



spin 



V,; 



(19) 



G=corasf 

This means that Q S pin{ G i\ is the effective spin Hamiltonian. The last term in Eq.flTSp is 
nothing but Dzialoshinskii- Moriya interaction term. It is non-zero only in relativistic case 
where E^ and can be, generally speaking, "non-parallel" and Gij ^ Gji for the crystals 
without inversion center. In the following we will not consider this term. 



C. Exchange interactions 



In the nonrelativistic case one can rewrite the spin Hamiltonian for small spin deviations 
near collinear magnetic structures in the following form 



n 



spin 



(20) 



'j 



where 



Jij — —Tr^L \^\G\^ a jG\i) 



(21) 



are the effective exchange parameters. This formula generalize the LSDA expressions of |7j 
to the case of correlated systems. 

The sum rule (Eq. |16D for the collinear magnetic configuration takes the following form 



"ii "it 2—d ij j j 



(22) 



Using Eq. ([22]) we obtain the following expression for the total exchange interaction of a 
given site with the all magnetic environment. 



vs/it v s ri v s ( rA n 



(23) 
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Spin wave spectrum in ferromagnets can be considered both directly from the exchange 
parameters or by the consideration of the energy of corresponding spiral structure (cf 0). 
In nonrelativistic case when the anisotropy is absent one has 

^ = ^ E J oj (l - cos qR,) = A [ J(0) - J(q)] (24) 

where M is the magnetic moment (in Bohr magnetons) per magnetic ion. Corresponding 
expressions can be easily written in k-space also. In the short notation q = (q, iu) it is easy 
to write the general expression for J(q) = J(q, 0): 

J (l) = -TT Y, Tr Lp s (k)GHk)i: s (k + q)G^k + q)} (25) 
7V fc k 

where N). is the total number of k-points. 

It should be noted that the expression for spin stiffness tensor D a/3 defined by the relation 

^ q = D a f3q a qp (26) 

(q — > 0) in terms of exchange parameters has to be exact since it is the consequence of phe- 
nomenological Landau- Lifshitz equations which are definitely correct in the long- wavelength 
limit (cf 0). One has from Eqs. (0,0) 

,0G* (k)8G^ (k)' 



n - 2 V/r¥- 2 Tr v ^Gt ( k ) 



dk 



(27) 

1V1 \ UKiX UKfj J 

where k is the quasimomentum and the summation is over the Brillouin zone. In cubic crys- 
tals D aj 3 = D5 a p. For arbitrary q, the expression of magnon spectrum in terms of J^- is valid 
only in the adiabatic approximation, i.e. provided that the magnon frequencies are small 
in comparison with characteristic electronic energies. Otherwise, collective magnetic exci- 
tations which are magnons cannot be separated from non-coherent particle-hole excitations 
(Stoner continuum) [|Hj and magnon frequencies ( PH) are not the exact poles of transverse 
magnetic susceptibility (which are even not real at large q). 

Now we have to consider the accuracy of expressions for Q2~H) themselves. Eqs. ( |T2"D 
and ( \L% are exact in LDA++ approach (i.e. with the only assumption about the local 
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self-energy). Hence, if one postulate the existence of effective spin Hamiltonian in the sense 
of Eq. (fl9|) , Eq. ([H]) is also exact. However, they do not have rigorous connection with the 
static transverse spin susceptibility. The latter is expressed in terms of the matrix 



6 2 Q 5V« 



without the restriction G = const. They differ from our exchange parameters by the terms 
containing From the diagrammatic point of view in the framework of DMFT Q they 
are nothing but the vertex corrections. We do not present the corresponding expressions 
since the benefits of the introducing of exchange parameters beyond adiabatic approximation 
which is equivalent to (]T3) is doubtful. In more rigorous consideration it is convenient to 
work directly with DMFT expressions for static and dynamic susceptibility [|]]. 

Thus one can see that, generally speaking,the exchange parameters differ from the ex- 
act response characteristics defined via static susceptibility since the latter contains vertex 
corrections. At the same time, our derivation of exchange parameters seems to be rigor- 
ous in the adiabatic approximation for spin dynamics when spin fluctuation frequency is 
much smaller than characteristic electron energy. The situation is similar to the case of 
electron-phonon interactions where according to the Migdal theorem vertex corrections are 
small in adiabatic parameter (ratio of characteristic phonon energy to electron one) [12 



The derivation of the exchange parameters from the variations of thermodynamic poten- 
tial, being approximate, can be useful nevertheless for the fast and accurate calculations of 
different magnetic systems. 

Note that in the LDA++ approach, as well as in LDA+U method and in contrast with 
usual LSDF, one can rotate separately spins of states with given orbital quantum numbers 
L, V . For example, for non-relativistic case one can obtain 

^spin ^ ] JiLjjL'&iL^jL' 

iL,jL> 

where 

JiL,jL' = —Tr^ {^iLG\ L ^ L iY^ L ,G^ L , . 
9 



are orbital dependent exchange parameters. 



D. Magnetic anisotropy 

Let us consider now the change of spin energy at the rotation of all the spins at the same 
angle. It is definitely zero in nonrelativistic case. In the presence of spin-orbit coupling, it 
is nothing but the energy of magnetic anisotropy. One can obtain from Eq.fll8"D 

tt ams = Tr wL £ { (G?- x £*) (G£ x Ef ) } = (28) 
Tt^l E {( GS ( k ) x SS ) (G S (-k) x S s )} 

k 

where the last equality is valid for ferromagnets with one magnetic atom per unit cell, k is 
the quasimomentum and the summation is over Brillouin zone. 

Finally, note that we use essentially three properties of LDA++ approach: (i) locality 
of self-energy (ii) spin-independence of bare Green function (i.e. spin-independence of bare 
LDA spectrum; all magnetic effects including Hartree-Fock ones are included in self- energy) 
and (iii) approximations for self energy have to be conserving, or "$- derivable" since only 
in that case analog of "local force theorem" (||) takes place. 



III. EXCHANGE INTERACTIONS IN FERROMAGNETIC IRON 

A. Computational technique 

As an example we calculate the magnetic properties of ferromagnetic iron using the most 
accurate method to take into account local correlations. For this purpose we use the local 



quantum Monte-Carlo approach B with the generalization to the multiband case [13 
We start from LDA+U Hamiltonian in the diagonal density approximation: 

H = E tim,fm' C trna C i'rn'cr + (29) 
{ima} 

^ , U mm rTli mo -Tli m i _<j + 



2 

vmm' a 
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2 ^ mm' J mm') n ima n im' a 
imj^m'cr 

where t LDA is effective single-particle Hamiltonian obtained from the non-magnetic LDA 
with the corrections for double counting of the average interactions among d-electrons; i is 
the site index and m is the orbital quantum numbers; a =|, { is the spin projection; c + ,c 
are the Fermi creation and annihilation operators (n = c + c). 

The screened Coulomb and exchange vertex for the d-electrons 

Umm> = < mm' I V s e ^ r (r - r') \mm! > (30) 
J mm' = < mm'\VZ(r - r')\m'm > 



scr\ 



are expressed via the effective Slater integrals and corresponds to the average U = 2.3 
eV and J = 0.9 eV ||. We use the minimal spd-b&sis in the LMTO-TB formalism and 
numerical orthogonalization for t LDA (k) matrix ||. The Matsubara frequencies summation 
in our calculations corresponds to the temperature of about T=850 K. 
Local Green-function matrix has the following form 

G(zu) = ^{iu + t LDA (k) - S(^)}- 1 (31) 

k 

Note that due to cubic crystal symmetry of ferromagnetic bcc-iron the local Green func- 
tion is diagonal both in the orbital and the spin indices and the bath Green function is 
defined as 

G° m (ico) = G m (ico) + E m (iu) (32) 

The local Green functions for the imaginary time interval [0, /3] with the mesh t\ = I At, 
I = 0, ...,L — 1, and At = (3/L, where (5 = ^ is calculated in the path-integral formalism 

11B 



G£ = ! J2det[0(s)]*G l i(s) (33) 

mm' 

here we redefined for simplicity m = {m,a},Z is the partition function and the so-called 
fermion-determinant det[0(s)] and the Green function for arbitrary set of the auxiliary 
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fields G(s) = O 1 (s) are obtained via the Dyson equation [13] for imaginary-time matrix 
(G m ( a )=G«[( S )): 



G m =[l-(G^-l)(e v '--l)]- 1 G 



where the effective fluctuation potential from the Ising fields s l mm , = ±1 is 



1, m < m' 
— 1, m > w! 

and the discrete Hubbard- Stratonovich parameters are \ mm < = arccosh[exp(iAr[/ mm /)] 
The main problem of the multiband QMC formalism is the large number of the auxiliary 
fields s l mm i. For each time slices I it is equals to M(2M — 1) whereM is the total number 
of the orbitals which is equal to 45 for d-states. We compute the sum over this auxiliary 
fields in Eq.|33] using important sampling QMC algorithm and performed a dozen of self- 



consistent iterations over the self-energy Eqs. (pT] . [32| . |33D . The number of QMC sweeps was of 
the order of 10 5 on the CRAY-T3e. The final G m (r) has very little statistical noise. We use 
maximum entropy method [15] for analytical continuations of the QMC Green functions to 
the real axis. Comparison of the total density of states with the results of LSDA calculations 
(Fig.0) shows a reasonable agreement for single-particle properties of not "highly correlated" 
ferromagnetic iron. 

Using the self-consistent values for H m (iu)) we calculate the exchange interactions (Eq.pB]) 
and spin-wave spectrum (EqJ24j) using the four-dimensional fast Fourier transform (FFT) 
method for (k, iui) space with the mesh 20 3 x 320. We compare the results for the 
exchange interactions with corresponding calculations for the LSDA method 0. 



12 



B. Computational results 



The spin-wave spectrum for ferromagnetic iron is presented in Fig.0 in comparison with 



the results of LSDA-exchange calculations and with different experimental data [16-18 



Note that for high-energy spin- waves the experimental data [18] has large error-bars due to 
Stoner damping (we show one experimental point with the uncertainties in the "g" space). 
On the other hand, the expression of magnon frequency in terms of exchange parameters 
itself becomes problematic in that region due to breakdown of adiabatic approximation, 
as it is discussed above. Therefore we do not believe that good agreement of LSDA data 
with the experimental ones for large wave vectors is too important since it may be a result 
of "mutual cancellation" of inaccuracies of adiabatic approximation and LSDA itself. At 
the same time, our lower-energy spin-waves spectrum (where the theory is reliable) agree 
better with the experiments then the result of LSDA calculations. Our LSDA spin-wave 



spectrum agree well with the results of frozen magnon calculations [19,20]. Note that in 



the LSDA scheme one could use the linear-response formalism J2TJ to calculate the spin- 
wave spectrum with the Stoner renormalizations, which should gives in principle the same 
spin-wave stiffness as our LSDA calculations. The corresponding exchange parameters and 
spin-waves stiffness (Eq.|26[) are presented in the Table I. The general trend in the distance 
dependence of exchange interactions in ferromagnetic iron is similar in both schemes, but 
relative strength of various interactions is quite different. Experimental value of spin- wave 
stiffness D=280 meV/A 2 ][17| agrees well with the theoretical LDA++ estimations. 



IV. CONCLUSIONS 

In conclusion, we present a general method for the investigation of magnetic interactions 
in the correlated electron systems. This scheme is not based on the perturbation theory in 
"[/" and could be applied for rare-earth systems where both the effect of the band struc- 
ture and the multiplet effects are crucial for a rich magnetic phase diagram. Our general 
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expressions are valid in relativistic case and can be used for the calculation of both ex- 
change and Dzialoshinskii- Moriya interactions, and magnetic anisotropy. An illustrative 
example of ferromagnetic iron shows that the correlation effects in exchange interactions 
may be noticeable even in such moderately correlated systems. For rare-earth metals and 
their compounds, colossal magnetoresistance materials or high-T c systems, this effect may 
be much more important. For example, the careful investigations of exchange interactions 



in MnO within the LSDA, LDA+U and optimized potential methods for MnO show 
the disagreement with experimental spin-wave spectrum (even for small q) , and indicate a 
possible role of correlation effects. 

As for the formalism itself, this work demonstrates an essential difference between spin 
density functional approach and LDA++ method. The latter deals with the thermodynamic 
potential as a functional of Green function rather than electron density. Nevertheless, there 
is a deep formal correspondence between two techniques (self-energy corresponds to the 
exchange- correlation potential, etc). In particular, an analog of local force theorem can 
be proved for LDA++ approach. It may be useful not only for the calculation of magnetic 
interactions but also for elastic stresses, in particular, pressure, and other physical properties. 
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TABLES 

TABLE I. Parameters of exchange interactions and spin- wave stiffness for ferromagnetic iron 
calculated with the LSDA and LDA+S methods. 



meV 


Jo 


Ji 


h 


J 3 


J 4 


h 


J 6 


D (meV/A 2 ) 


LSDA 


166.1 


16.48 


8.07 


0.25 


-1.03 


-0.31 


0.26 


320 


LDA+S 


115.8 


13.31 


2.5 


0.73 


-0.38 


-0.83 


0.01 


260 
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FIGURES 
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Energy ( eV ) 



FIG. 1. Total spin-polarized density of states for ferromagnetic iron in the LSDA and 
LDA+QMC approximations. 
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FIG. 2. The spin-wave spectrum for ferromagnetic iron in the LSDA and LDA+S approxima- 
tions compared with different experiments (circles [16], squares [17], and diamonds [18]) (a); The 
corresponding spin- wave spectrum from LDA+S scheme in the (110) plane (b). 
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